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Aims. We investigate Galactic bulge planetary nebulae without emission-line central stars for which peculiar infrared spectra have been 
obtained with the Spitzer Space Telescope, including the simultaneous signs of oxygen and carbon based dust. Three separate sub- 
groups can be defined characterized by the different chemical composition of the dust and the presence of crystalline and amorphous 
silicates. 

Methods. We use literature data to analyze the different nebular properties and deduce both the evolutionary status and the origin 
of these three groups. In particular, we check whether there are signs of evolutionary links between dual-dust chemistry planetary 
nebulae without detected emission-line central stars and those with emission-line stars. 

Results. Our primary finding is that the classification based on the dust properties is reflected in the more general properties of these 
planetary nebulae. However, some observed properties are difficult to relate to the common view of planetary nebulae. In particular, it 
is challenging to interpret the peculiar gas chemical composition of many analyzed objects in the standard picture of the evolution of 
planetary nebulae progenitors. We confirm that the dual-dust chemistry phenomenon is not limited to planetary nebulae with emission- 
I . line central stars. 

Conclusions. Our results clearly indicate that there is no unique road to the formation of planetary nebulae even in a homogeneous 
environment such as the Galactic bulge. The evolution of a single asymptotic giant branch star may lead to the formation of different 
types of planetary nebulae. In addition, the evolution in a close binary system should sometimes also be considered. 
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(Nj ■ 1- Introduction rich). Features at 3.3, 6.2, "7.7", 8.6, and 11.3/an attributed to 

\/-\ polycyclic aromatic hydrocarbons (PAHs) are common in C-rich 

J* > ; After the completion of hydrogen and helium burning in their pNe { Garcfa . Lario et aL 1999i and references therein) while 

^ , cores, low- to intermediate-mass stars (0.8 < M < 8 Mp) evolve g features attributed t0 crystanine silicates (e.g., those cen- 

. ■ towards the asymptotic giant branch (AGB; e.g., | Herwig || 2005|) t ered on 23.5, 27.5 a nd 33.8 fim) are usually found in O-rich PNe 

Cn ; and then pass through the planetary nebula (PN, plural PNe) ^ Molster et al 2002|) 

— phase before ending their lives as white dwarfs. At the tip of 1 c „ , .. ,. , 1 XT , , , , , , , 

, . , , . , , A few Galactic disk PNe exhibited a remarkable dual-dust 
the AGB phase, these stars experience strong mass loss that 
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... . , . , , .. ,, ,. , (C-rich and O-rich) chemistry showing both PAH and crys- 

efficiently enriches the surrounding interstellar medium with features in ISO snertral Warers eta l 1998ab- 

>. j. . j A ~ . - . nr , . . Ldiiine siiicdLe ledmies in lovj specua ^waieis ei ai. iwoa.D, 

huge amounts of gas and dust. Stars leave the AGB when the tttt — nr: nnn « ml , m. c r . . T~- * 

. _h , , , , j, , , , Cohen etal. 1999, 2002). The fact that this was an infrequent 

strong mass loss stops and then the future central star (CS, plural — 7" — T . , ,. , . 

JC* . „„ . -j, , . ■ , , rr -i phenomenon may be due to the instruments used, which in 

rN , CSs) rapidly evolves towards hotter effective temperatures in the r , , , ^ ,,. ... 

\-k TT : t, „ j- rr, , .... „ . many cases may have been unable to detect crystalline silicates. 

rrt ■ Hertzsprung-Russell diagram. Thus, when the ionization of the „ J , v ' „ . „ „ . ,„ rrr, — n 

tvJ . . f i , r,»r • r i tt • For example, the Spitzer Space Telescope (Spitzer, Werner et al. 

elected gas takes place, a new PN is formed. However, in most l»„a , \ . , r . ... .,. . . ,. . 777^ — 

J ^ . . , j- ■ ■ , • ,, ■ 2004), detected crystalline silicates in 10 post-AGB sources 

cases the total amount of ionized gas is very small compared to i„ ■' 1 -i ^ , • , 

. ... . A • ' , , ,. (Cernaone et al. 2009) while after completion of the ISO mis- 

the total mass previously elected. An important traction of this . , — z ~r . ^ — I — TIRnm \ 

. , . : f. J , . j. 1 . , , sion only 2 such sources were known (Szczerba et al. 2003). 
material remains neutral in the form of dust grains, molecules, _ . , , . „ XT ,. — 7^ — ~~~ — — : — , 
or atoms, which can be easily detected in the infrared domain. , The mixed-chemistry PNe discovered by ISO pertain to the 
Thanks to the analysis of the features observed by the Infrared C ^ 1 SS J °L?^ e r ; t x S r Wlth u C " u nch Wolf-Rayet type nuclei (the so- 
Space Observatory (ISO) in the spectra of PNe, it was possible c u alled [WR] PNe), which usually show a lack of hydrogen in 
to confirm the presence of large amounts of dust grains around their atmospheres. These atmospheres are instead mostly com- 
PNe as well as the dominant dust chemistry (C-rich versus O- P osed , of heliu m carbon, and o xygen and the CSs show intense 
J mass-loss (e.g.. lCrowtherll2008l) . 

Send offprint requests to: S.K. Gorny The evolution of an AGB star with a stellar or substellar 

* Based on observations made with the Spitzer Space Telescope, companion that undergoes the common-envelope phase is an- 

which is operated by the Jet Propulsion Laboratory, California Institute other possible way of creating a PN. Some authors argue that a 

of Technology, under NASA contract 1407. companion object is often mandatory for a planetary nebula to 
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be created (see in Ide Marcoll2009t and references therein). We 
noted that some hypotheses compiled to explain the simultane- 
ous presence of carbon and oxygen dust also require a binary 
system. 

The Galactic bulge is characterized b y an old population 
of mostly low-mass st ars (IZoccali et al.l {2003), but see also 
lUttenthaler et al.l J2007) and references therein). It is also known 
that the abundances of PNe in the Galactic bulge (GBPNe) 
differ from those located in the Galactic disk as they have 
higher metallicities and lower C/O ratios (e.g., I Wang & Liu! 
l2007h . The differences in metallicity seems to play a domi- 
nant role in the chemical evolution o f low- t o intermediate-mass 
stars (e .g.. iGarcfa-Hernandez et al.1 (l2007h : IStanghellini et al.l 
d2007l) : IChiappini et all d2009h )~Thus. studding the Galactic 
bulge enables us to investigate the stellar evolution of low- and 
intermediate-mass stars in higher metallicity environments and 
at the same time an insight into the chemical evolution and for- 
ma tion of our Milky Way. 

iGutenkunst et alJ (12008) analyzed Spitzer spectra ac quired 
using the Infrared Spectrograph (IRS. iHouck et alj|2004h of 11 
PNe towards the Galactic bulge and inferred dual-dust chemistry 
in 6 of them. They suggested that the high percentage of dual- 
dust chemistry sources is unsurprising because the fraction of 
[WR] PNe is significantly higher in the bulge than in the Galactic 
disk. However, as we checked, only one of their dual-dust chem- 
istry sources have the wind characteristics of [WR] type CS 
and the higher proporti on of genuine [WR] PNe in the Galactic 
bulge is not confirmed (Gorn vet alJ2009l) . lPerea-Calder6n et al.l 
(2009) found that dual-dust chemistry is truly widespread among 
GBPNe. They analyzed a larger sample of 26 GBPNe observed 
with Spitzer/IRS among wh ich 21 exhibit dual-dust chemistry. 
iPerea-CaTd eron et al. (2009) observations shown that the simul- 
taneous presence of oxygen and carbon-rich dust features in 
the infrared spectra of [WR] PNe is not restricted to objects 
with late/cool [WC] class stars. In addition, dual-dust chemistry 
was seen in all observed PNe with WEL stars {weak emission- 
line stars, iTvlenda et al.lll993l) as well as members of recently 
discovered VL group (low io nization PNe arou nd stars with 
very late fWC 1 l]-like spect ra, iGornv et alj|2009l) . Surprisingly, 
IPerea-Cal deron et al. (2009) found dual-dust chemistry also in 
some PNe without detected emission-line CSs. 

Another interesting pr operty of the PNe observed by 
IPerea-Cald eron et al. (2009) was the amorphous silicate emis- 
sion at lOyum, which was detected in four dual-dust chem- 
istry GBPNe and in most of the O-rich PNe that they ob- 
served. Note that before Spit zer there was known only one such 
PN, namely SwSt 1 (e.g., ISzczerba etail 1200 lb . belonging to 
the fWR] PNe. In contrast, the 10/mi feature objects found by 
Perea-Calderon et al. (2009), neither belong to this group nor ex- 
hibit stellar emission lines. 

In this work, we analyze the multiple properties of PNe with- 
out emission-line CSs that are found to have peculiar infrared 
spectra acquired with Spitzer/IRS. The paper is organized as fol- 
lows. In Sect. 2 and 3 we describe our working sample and the 
main properties of their infrared spectra, respectively. We ana- 
lyze the nebular properties and evolutionary status of these PNe 
in Sect. 4. Finally, in Sect. 5 we discuss the results obtained. Our 
concluding remarks are given in Sect. 6. 

2. Sample selection 

We analyze the various properties of PNe without emission-line 
CSs observed with Spitzer/IRS and exhibiting signs of dual- 
dust chemistry and/or lOyum emission band of amorphous sil- 



Table 1. List of analyzed PNe. The Galactic bulge (b) or disk 
(d) association of the object is marked in Col. 3. The reference 
to original observers is given in Col. 4. Objects are divided into 
three groups according to their infrared properties as marked in 
Cols. 5-7 (see description in Sect. 3). 
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Ref: G - Gutenkunst et al. (2008); P - Perea-Calderon et al. (2009) 



icates. These PNe are listed in Table 1 along with a short 
descriptio n of their infrared fea t ures. Most of them were ob- 
served bv IPerea-Calderon et al.l (120091) . while two come from 
Gutenkunst et alJ (120081) . 

All of the sel e cted PNe were included in the analy- 
ses of IGornv et"al1 d2009l) but none were attributed to either 
WEL, VL, or [WR]-type groups of GBPNe. The absence of 
emission-line CSs was perceived in these PNe by studying 
high qualit y optical spect ra acqu ir ed at 2 and 4-me t er tele - 
scopes bv ICuisinier et al.l d2000l). lEscudero & Costal d200ll) . 



lEscudero et alJ d2004t). IGornv et al.l (120041) . and new observa- 
Gornv et all hooW/ 



tions in 



or b y checking the list of the ob- 
served lines in I Wang & Liul ( 20071). The s earche d stellar emis- 
sion lines were the same as in Go rny et all d2004l) . In none of the 



thirteen PNe presented in Table 1 did we notice any characteris- 
tics of emission-line CSs during additional direct re-inspection 
of the spectra. 

The objects for which the CS was of unknown type are 
not considered until their spectral classification is established. 
The PNe with dual-dust chemistry but unknown spectral type 
of the CS are H 2-20 a nd M 2-5 observed by Guten kunst et al. 
(200813 and H 1-62 bv IPerea-Calderon et alJ d2009l) . The spec- 
tral type of the CS is also unknown for three O-rich PNe with 
the detected lOji m emission feature, i.e., He 3 -1357, Cn 1-3 and 
IC 4732 listed bv lPerea-Calderon et alJ d2009h . 

Eleven of the PNe listed in Table 1 have a high probability of 
physically belonging to the Milky Way bulge because they sat- 
isfy the standard criteria (Stasinska & Tylenda 1994), namely: a) 
they are located within 10 degrees of the center of the Galaxy, b) 
have diameters smaller than 20", and c) the known radio fluxes 
at 5GHz are lower than 100 mJy. The two remaining objects are 
either too bright in the radio domain (Hb 6) or marginally out- 
side the 10 degree angular distance (He 2-260) and therefore 
regarded as possible Galactic disk members. 



1 The remaining dual-dust chemistry PNe on their list belong to [WR] 
type (M 2-3 1 ) and WEL PNe (H 2- 1 1 ). 
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The 13 objects collected in Table 1 are compared in this pa- 
per with the remai nder of the total of 180 GBPNe analyzed in 
iGornv et"ail d2009t) . Among this large reference group, there are 
119 PNe without emission-line CSs (which we refer to as normal 
PNe throughout this paper), 25 WEL, 14 VL, and 9 [WR] PNe. 
Some of the PNe from the latter groups were also observed with 
Spitzer/IRS by either the Perea-Calderon or Gutenkunst teams. 
This includes 5 WEL, 4 [WR], and 3 VL PNe. 

In this paper, we compare the objects investigated here to 
both normal PNe and PNe with emission-line CSs since it can 
not be excluded that the latter, in particular [WR] PNeQ, have at 
times weaker stellar winds and emission lines from CSs are not 
observed. It has to be checked if this is not the case of some of 
the PNe from Table 1. On the other hand, the stellar emission 
lines of WEL and VL PNe are not very strong and may escape 
unnoticed in spectra of lower signal to noise. Therefore, to con- 
firm that indeed dual-dust chemistry phenomenon also occurs in 
PNe without emission-line CSs we compare their properties with 
those of [ WR] , WEL, and VL PNe. 

Among the largest group of 1 19 normal bulge PNe that we 
use as a reference sample only four have been observed with 
Spitzer/IRS. This is a clear result of the s election effect since 
the tar gets for observations were chosen bv lPerea-Calderon et al.l 
(2009rbased on high quality IRAS fluxes at 12, 25, and 60 pm. 
In particular, the requirement of reliable measurements at 12/im 
is fulfilled by only 15 normal PNe of 64 with IRAS data avail- 
able. This means that no prominent emission features around 
12yum should be expected in normal PNe. Even in the case of 
the four objects observed by Spitzer, only crystalline silicates are 
present in their spectra and both PAH and amorphous silicates 
seem to be absenfl Unfortunately, the existing optical spectra of 
these four PNe are at the same time of low quality implying e.g., 
that the nature of their CSs is uncertain and many parameters 
cannot be reliably determined. 

We add to the discussion the analysis of two PNe from the 
Galactic disk population. One is the [WR]-type SwSt 1 men- 
tioned already in the introduction and the other is IC4776 
that has a WEL CS and was observed by|P erea-Calderon et al. 
(2009). Both these objects exhibit the dual-dust chemistry sig- 
natures typical of PNe with emission-line CSs but additionally 
also the IQ/im amorphous silicate feature. 

3. Infrared spectra 

The infrared spectra of the GBPNe withou t emission-Une CSs 
selecte d for this work were acquired by IPerea-Calderon et"aT] 
(120091) and iGutenkunst et alj d2008l) . The description of the 
Spitzer/IRS data reduction process is not repeated here and we 
refer the reader to these two references for details. 

After a short inspection of the Spitzer/IRS spectra, we find 
that among PNe without emission-line CSs, there are three 
clearly different types of objects. The examples are presented 
in Fig.Q] It can be seen that H 1-16 and Th 3-4 PNe have mixed 
chemistry where both PAHs and crystalline silicates features are 
detected. For Th 3-4 we also detect the amorphous silicate emis- 
sion feature at 10 fim. On the other hand, M2-23 is clearly O- 
rich (there is a lack of carbon-based dust features such as PAHs), 



2 In the Galactic bulge, only [WR] PNe with CSs of intermediate 
[WC] spectral classes are known, in co ntrast to the situation in the 
Galactic disk (see e.g. . [Go rny et al. 20011) 

3 Among these PNe, object M3-13 has exceptional Spitzer spectra 
with unide ntified broad emission feat ures at shorter wavelengths, see 
Fig.B4 of lPerea-Calderon et alj d2009h 



exhibiting the amorphous silicate feature at 10 pm and weak 
crystalline silicate features at longer wavelengths. Finally, we 
note that M 2-23 contains a circumstellar dust envelope, which 
is hotter than in H 1-16 and Th 3-4 as can be inferred from the 
shape of continuum emission in the spectrum. 

Based on the examples from Fig.Q] we can describe the spec- 
tra of all PNe analyzed in this work and divide them into differ- 
ent subgroups (Table 1). Each group has some similarity with 
one of the other groups and clear dissimilarity with another. The 
first group of 5 PNe listed at the top of Table 1 is character- 
ized by the simultaneous presence of both carbon-based dust 
(PAH features at 6.2, "7.7", 8.6, and 11.3 pm) and oxygen-based 
dust (crystalline silicate features at 23.5, 27.5, and 33.8/mi). 
We refer to this group as DC cr (dual-dust chemistry with sili- 
cates only in crystalline form) subsample throughout the rest of 
this paper. The second group in Table 1 has 4 members and is 
also characterized by the simultaneous presence of oxygen and 
carbon dust but at the same time there are also signs of amor- 
phous silicate features at about 10/um. We therefore refer to this 
group as DC a+cr (dual-dust chemistry with silicates in amor- 
phous and crystalline form). Finally, 4 PNe in Table 1 exhibit 
only oxygen dust features, which however include the uncom- 
mon lOyum amorphous silicate feature. This last group is called 
OC fl+c , (oxygen-dust chemistry with amorphous and crystalline 
forms) subsample. 

Detailed comparison and analysis of dust features observed 
with Spitzer/IRS in GBPNe will be presented in the forthcoming 
paper (Szczerba et al. in prep.) 



4. Nebular properties and evolutionary status 

We analyze different nebular properties that can shed the light 
on the three groups of bulge PNe and allow us to infer their 
evolutionary status and origin. Since some of them (in partic- 
ular DC tT ) have infrared spectra resembling those of PNe with 
emission-line CSs, we also check whether any other properties 
would associate them with such objects. 

In Table 2 we report the basic nebular parameter s of the in- 
vestig ated PNe. Most of them were adopted from IGornv et al.l 
(200 91), recomputed using the methods and formulae from 
Stasiriska & Szczerbal ([1999), or adopted from other literature 
sources. 

4.1. Bulge location 

As mentioned in Sect. 2, most of the PNe analyzed here are not 
only simply observed towards the center of the Milky Way but 
most probably physically pertain to the Galactic bulge. In partic- 
ular, they are located less than 10° from the center of the Galaxy. 
Figure|2]presents the radial velocities V/„. of bulge PNe corrected 
for solar motiorQ as a function of Galactic longitude coordinate 
/. Although the samples of DC„-, DC fl+cr , and OC fl+e) - are not nu- 
merous, we can state that members of each of them are found at 
different longitudes within the bulge and no grouping at special 
locations can be distinguished 



4 The rad i al vel ocities used in this paper have been taken from 
iDurand et all {l9?8) and corrected to the local standard of rest using 
the formulae of Beaulieu et al. ( 200G) 

5 Recently Gorny et al. ( 2009) noticed a difference in locations be- 
tween bulge [WR], VL, and WEL PNe with the last two groups typi- 
cally being located at longitudes less than 4° 5 from the center whereas 
bulge [WR] PNe are often found at greater longitudes. Obviously, the 
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Fig. 1. Examples of Spitzer/IRS spectra of PNe with three different types of dust composition: H 1-16 of DC cr type at the top, 
Th3-4 of DC, +cr in the middle, and M 2-23 of OC fl+(T type at the bottom. Dotted vertical lines indicate the positions of crystalline 
silicate emission features at 23.5, 27.5, and 33.8 fim (left panels) and typical PAH features at 6.2, "7.7", 8.6, and 11.3 /mi (right 
panels). The lO/vm amorphous silicate feature is seen in the middle and bottom spectra. 



As can be noted in Fig.|2l the velocities of Galactic bulge 
DC cr , DC a+cr , and OC a+cr objects (large filled symbols) are usu- 
ally large or very large and therefore their kinematic properties 
are typical of PNe physically related to the bulge system (com- 
pare with Fig. 13 of IGornv et alJ ([2004)). On the other hand, the 
two objects already assigned to Galactic disk (Hb 6 and He 2- 
260) have small V/ sr velocities that are characteristic of that PNe 
population. 

For the other Galactic coordinate, the latitude b, the bulge 
PNe analyzed here show nothing peculiar in their distribution. 
They are found at b smaller than about 5°, but, as for all the other 
known PNe, avoid latitudes below 2°. This is simply because 
of the interstellar dust prohibiting their discovery. However, in 
Fig.[3]we plot the bulge PNe in b coordinates versus logarithmic 
extinction parameter at Hfi (C opt ) and some display a very in- 
teresting property: the members of DC fl+c) (in particular Th 3-4 
and H 1 -40) have much greater extinction than expected of typi- 
cal PNe located at their latitudes off the Galactic plane. 

The extinction plotted in Fig. [3] (see also Col. 2 of Table 2) 
were derived from the ratios of hydrogen lines in optical 
spectra. This was accomplished mainly by comparing the ob- 
served Balmer HafH/3 ratio with its theoretically expected value, 
although sometim es Ha/Hy had to be used (see details in 
IGornv et"al]|2009t) . We checked that when extinction can be de- 

4! 5 longitude does not seem to be important for any of the subsamples 
investigated here. 



rived from both of these ratios the agreement is usually very 
good. 

The value of PNe extinction may be above average when 
the matter along the line of sight has a different blocking prop- 
erties or if the object is simply located behind a larger amount 
of interstellar dust. This is the most straightforward explanation, 
although the chances of finding all four DC a+cr with these con- 
ditions can be evaluated at only a few percent. Therefore at least 
partly the excessive extinction towards DC fl+e) . could also be at- 
tributed to some specific properties of their own dust or the dust 
nearby around them. In that way, DC a+c) would be interesting 
candidates for studying the internal extinction of PNe. 

4.2. Diameters and densities 

In Fig. [4] the reddening-corrected fluxes of the hydrogen H/3 line 
of GBPNe are presented as a function of apparent angular di- 
ameters (Col. 6 in Table 2). As a reference we overp lotted the 
theoretical tracks for PNe at a distance of 8.5 kpc with CSs of 
three different masses 0. 57, 0.60, a nd 0.6 4M(T), evolving accord- 
ing to the predictions of iBlockerl d 19951) and assuming they ra- 
diate like a black body. The surrounding nebula were described 
by a simple model of uniformly filled sphere with total gas mass 
of O.2M0, filling factor e=0.75, and expansion at the constant 
velocity of 20km/s. As can be seen in this figure, the DC cr and 
DC a+cr PNe occupy a rather restricted region of the plot. Their 
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Table 2. Observational data for analyzed planetary nebulae. The reported quantities are: logarithmic extinction at H/3 (Col. 2); 
ionization parameters ++ /(0 + +0 ++ ) and He ++ /(He + +He ++ ) (Col. 3 and 4); electron density from [S II] A67 17/6731 line ratio 
(Col. 5); nebular diameter (Col. 6); nebular expansion velocity (Col. 7); kinematic age (Col. 8); H/3 surface brightness (Col. 9); 
stellar temperature from Zanstra hydrogen method (Col. 10); mass of ionized nebular gas (Col. 11) and mass of dust (Col. 12). 
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Fig. 2. Radial velocity versus Galactic longitude coordinate of 
the investigated PNe. The big filled symbols mark: DC„. as 
red hexagons, DC fl+e/ - brown diamonds, and OC fl+c , as violet 
reversed triangles. Galactic bulge PNe with emission-line CSs 
are marked with smaller open symbols: [WR] as pink trian- 
gles, WEL as blue squares, and VL as green circles. Normal 
Galactic bulge PNe are presented with small dark grey dots. 
Rotation curves for PNe at 1, 4, and 6 kpc circular orbits as- 
suming Galactocentric rotation velocities of 220 km/s are plotted 
with solid, dashed, and dotted lines respectively. 



H/3 fluxes are typically 2 times lower than those of the majority 
of [WR] PNe in the bulge, comparable to the brightest of WEL 
PNe and usually brighter than VL PNe. The diameters of DC„ 
and DC a+cr PNe range from about l"to 5", suggesting that they 



Fig. 3. Locations of investigated PNe in the Galactic latitude 
coordinate versus logarithmic extinction at H/3 parameter C opt . 
The meaning of the symbols is the same as in Fig. [2] 



are young or the nebulae are expanding slowly. However, expan- 
sion velocities are known for most of these objects (see Col. 7 
of Table 2) and seem norma l with a typical value o f 20 km/s. 
Using the data collected in iGesicki & Ziilstral (120071) . one can 
check that for the remaining GBPNe, normal, [WR], WEL and 
VL PNe combined (37 entries in their Table 3), the expansion 
velocities have a rather flat distribution from 10 to 32 km/s with 
a median value of V exp -22 km/s. 

The diameters of OC fl+CT PNe are typically only 1", which 
are the smallest among known GBPNe. However there are in- 
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Fig. 4. The relation between apparent diameter and reddening- 
corrected Hy6 flux for Galactic bulge PNe. The lines present 
model calculations for CSs of O.57M0 (dashed lines), O.6OM0 
(dotted), and 0.64Mq (solid), adopting the total nebular mass of 
0.20Mq. The meaning of the symbols is the same as in Fig. [2] 
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Fig. 5. Distributions of electron densities for the different groups 
of Galactic bulge PNe. For the normal PNe (bottom) the median 
value along with the 25 and 75 percentiles are marked with three 
short vertical lines above the histogram. Total numbers of objects 
included are shown in the left-hand parts of the panels below 
sample names. 



dications, that their expansion velocities may be smaller than 
average as they range from lOkm/s to only 17km/s. 

Assuming a simple model of constant nebular expansion that 
starts when the object leaves the AGB and knowing the angular 
diameter, distance, and V exp , it is possible to derive the kine- 
matic ages tki„ of PNe. They are listed in Col. 7 of Table 2 and 
range from 500 to 3800 years. For comparison , the median age of 
GBPNe is 2900 years using data collected in lGesicki & Ziilstra 
(I2007M This means that the objects analyzed here may not be 
extremely young but clearly belong to the younger PNe in the 
Galactic bulge. This is unsurprising taking into account the pos- 
sible selection effects when choosing objects for Spitzer obser- 
vations, as described in Sect. 2. 

In Fig. [5] we present the distributions of electron densities 
derived from [SII] 67 17/31 A line ratio for the different groups 
of GBPNe. The values of n e for individual objects can be found 
in Col. 5 of Table 20. Electron density is another independent 
parameter related to the age of the object. Surprisingly, as can 
be seen in Fig. [5] objects analyzed in this paper are among the 
densest of PNe in the Galactic bulge region. In particular, OC a+cr 
and DC U+CT PNe, i.e., the objects with the 10/im emission feature 
are much denser than the median log n e =3. 14 of normal GBPNe. 
The fact that they are also denser than the 75 percentile value of 
normal PNe implies that the distributions are truly different. The 
OC fl+cr and DC fl+ff are also considerably denser than the PNe 
with emission-line CSs. 



6 This median age of GBPNe may be underestima ted. There is a pos- 
sibility of important selection effects in the sample of Gesic ki & Ziilstral 
(2007) that seems strongly biased towards PNe with small angular di- 
ameters as the median 6 is 3'.'5 in their sample compared with 5'.'4 of 
our normal PNe. 

7 Table 4 is available online. 



4.3. Evolutionary status 

The evolutionary state of the central star of the planetary nebula 
is best described by its temperature. In Col. 10 of Table 2, we 
present the stellar temperatures estimated with the Zanstra hy- 
drogen method. Unfortunately, the data necessary to derive tem- 
peratures are not available for many DC cr and DC a+cr objects 
and it is then difficult to compare them with the temperatures of 
the other GBPNe. 

However, the most characteristic features of the planetary 
nebula spectrum are emission lines of many different atoms 
and their ions at different levels of excitation. As the tem- 
perature of the CS increases, the lines of increasingly higher 
ionization stages become observable. We take advantage of 
this property to at least qualitatively investigate the evolution- 
ary stage of all the objects, including those with unknown CS 
temperature. For this purpose we divide them into three ex- 
citation classes: i) those characterized by nebulae in which 
most of the oxygen atoms remain in the form of + ions 
i.e., O ++ /(O + +O ++ )<0.30 (for the coolest CSs); ii) those with 
most of their nebular oxygen in the form of ++ but no or 
very little He ++ ions present i.e., O ++ /(O + +O ++ )>0.30 and 
He ++ /(He + +He ++ )<0.03 (for the intermediate temperature CSs); 
iii) the remaining objects with considerable amounts of helium 
He ++ ions i.e., He ++ /(He + +He ++ )>0.03 (for the hottest CSs). 
Using the three excitation classes so defined, we checked sepa- 
rately for each group of PNe discussed in this work the distribu- 
tions of Hfi surface brightness. The parameters ++ /(0 + +0 ++ ) 
and He ++ /(He + +He ++ ) for individual DC„-, DC a+cr , and OC fl+ „. 
PNe are given in Cols. 4 and 5 of Table 2 and the values of Sh/i 
in Col. 9 of that Table. 

The Sh/3 parameter constitutes a good measure of the evo- 
lutionary advancement of the nebula and changes by a few or- 



Gorny et al.: New PNe with peculiar dust chemistry 



7 



ders of magnitude between the formation of the observable neb- 
ula and the moment it disperses into the interstellar medium. In 
Fig.|6l we present the histograms of Sh/3 for PNe with the coolest, 
intermediate, and hottest CSs in the left, middle, and right panels 
respectively. Above each histogram, the average nebular electron 
densities of the PNe that comprise the bars are overplotted. 

In Fig. [71 we present for comparison the theoretical predic- 
tions for the histograms of Sh/3- We used the same models as al- 
ready presented in Sect. 4.2 and Fig. [4] We considered PNe with 
CSs of three different masses, 0.57 0.60, and 0.64 M(T) evolv- 
ing according to the calculations of iBlockerl (Il995l) . The same 
simple model of surrounding nebula were assumed of uniformly 
filled sphere with total gas mass of O.2M0, filling factor e=0.75, 
and expanding at the constant velocity of 20km/s. To qualita- 
tively reconstruct the ionization conditions applied above to the 
real observed PNe, we divided the life of the model objects into 
three periods with CS temperature below 35 kK, between 35 kK 
and 70 kK, and finally above 70 kK (left, middle, and right panels 
in Fig. [7] respectively). The height of each bar is directly propor- 
tional to the period the object is characterized by the given Sfjp- 
It is normalized by the total time nebula should be bright enough 
to be detected (log F(H8)>-12.6) but not too extended and dis- 
persed (0<14arcsecfl The nebular electron densities expected 
from model calculations are also overplotted for the Sh/3 starting 
and ending limit of each histogram bar. 

Analyzing the distributions of DC„- and DC a+cr in Fig. [6] we 
note that PNe of both groups seem to be at a similar stage of neb- 
ular evolution. Roughly half of their population seems to be as- 
sociated with CSs of intermediate temperatures, while the other 
half has CSs that are already very hot. There are no DC cr nor 
DC„ +cr PNe with cool CSs but the apparent deficit of objects 
with cool CSs seems to be a general property of GBPNe (with a 
clear exception of the VL PNe). 

The DC„ and DC a+cr with stellar temperatures in the range 
35kK to 70kK are well reproduced by models with intermediate 
mass CSs (compare with predictions for 0.60Mq star in Fig.[7j). 
It is also important to note that their derived Sh/3 have values typ- 
ical of other GBPNe, including objects with emission-line CSs. 
Therefore, these DC cr and DC fl+c , PNe do not seem to be less 
evolutionary advanced than [WR] or WEL PNe and thus cannot 
be their predecessor^. 

A different situation can be seen in the right panel of Fig.|6l 
since the Sh/3 of DC cr and DC fl+cr are clearly larger than those of 
both normal PNe and of the few GBPNe with hot emission-line 
CSs. By comparing with the models in Fig. [7] it is apparent that 
DC C) - and DC a+cr PNe with CS temperatures above 70 kK origi- 
nate in more massive AGB progenitors than the group discussed 
above. Analyzing the model distributions, it can also be deduced 
that they should have a considerable number of unrecognized an- 
cestors among normal PNe with hot CSs. The evolutionary link 
to some PNe with emission-line CSs also cannot be excludec0- 

In summary, it seems that DC„- and DC a+cr PNe can be found 
among PNe with both massive and intermediate-mass CSs. As 
both groups are characterized by similar range of temperatures 
and Sh/3 values but have clearly different electron densities (see 
Fig. [5]), their evolutionary status cannot be identical. The de- 



8 Compare with the actual values for bulge PNe in Fig.[4] 

9 The existence of such predecessors with undiscovered or not 
yet active strong stellar winds could be one of the explanations why 
Galactic bulge [WR] PNe have only intermediate [WC] class CSs. 

10 The analysis is hampered by the Blocker ( 1995) tracks being ap- 
plicable directly only to PNe with H-burning nuclei, i.e., not to [WR] 
type CSs. 



tails may be derived by future complete modeling. However, the 
higher electron densities of DC n+c/ - PNe could be explained if 
the distribution of gas was more clumpy in these PNe. 

The Sh/3 distributions of OC a+ „- PNe can be compared with 
the separate theoretical predictions in Fig. I1;[PI Since the OC a+cr 
PNe evolve more slowly, the expansion velocity of 12km/s was 
assumed in these models (see in Sect. 4.4 and Table 2). In addi- 
tion, a total nebular mass of 0.1 OMq was adopted so the mod- 
els fit the small diameters of these PNe and do not exceed the 
measured electron densities. A qualitative agreement can be ob- 
served for OC a+cr PNe with models of intermediate-mass CSs 
while models assuming lower-mass progenitors can be rejected. 
The only exception is He 2-260, an object with the coolest CS 
and a nebula with the fastest expansion among OC a+cr PNe that 
is probably located in the Galactic disk. 

4.4. Masses of nebular gas and dust 

Since the objects that we analyze were selected based on the ap- 
pearance of their infrared spectra, it is important to learn about 
their dust content. The first important property is the dust-to-gas 
mass ratio md/m g . It can be easily estimated from observables by 
adopting a sim ple model of the dusty nebul a. We apply here the 
method from IStasinska & Szczerbal (Q~999) and recompute the 
paramet ers by adopting the new data on electron densities de- 
rived bv iGornv et alJ d2009l) . In addition, dealing with Galactic 
bulge objects we can adopt the distance of 8.5 kpc and derive ab- 
solute values of nebular ionized gas and dust masses. Our results 
are presented in Fig.[8]and in Cols. 1 1 and 12 of Table 2. For the 
Galactic disk PNe Hb 6 and He 2-260, the dust mass was calcu- 
lated from the derived md/m g ratio by assuming m g =0.2MQ. 

As can be seen in Fig. [8] the dust content of DC„ and DC a+cr 
PNe is usually similar to the typical value for normal GBPNe 
(median log;7?</=-3.30, represented by horizontal short-dashed 
line in the plot). It can be recalled here that DC fl+ „- can be 
suspected to have a considerable internal extinction (Sect. 4.1). 
However, since the dust content in these PNe is not exception- 
ally greater than in other objects, it could not explain excessive 
internal extinction. The reason for it would have to be in some 
specific properties of their dust that allow it to block more radi- 
ation than in other PNe. 

As far as the dust-to-gas ratio of the PNe presented in Fig. [8] 
is concerned, the DC a+cr represent the group with the highest 
ratio (see their location above the dotted line representing the 
median md/m g =2Al xl0~ 3 relation of normal PNe). The DC ff 
are much closer to this line meaning their md/m g ratios have nor- 
mal values. However, we recall that if some part of the nebula 
is still not ionized the value of ntd/m g that we derive may be 
overestimated. This is e.g., most probably t rue for the VL PN e 
that are in an extremely low ionization state (IGornv et al.l l2009). 
For DC a+cr PNe, the derived ionized gas masses are all below 
the median mass of normal PNe in Fig.|8](;« s; =0.18 Mq, marked 
with vertical long-dashed line). This indicates that the DC a+cr 
PNe may also be only partially ionized, if their total nebular 
masses are similar to those of normal PNe. 

As can be seen in Fig. [8] the mass of dust in individual 
OC fl+cr PNe may differ by more than an order of magnitude. For 
two of them M2-23 and H 1-32, we derived very small values 
of mj below 10~ 4 Mq. On the other hand, except for M2-23, 
the OC fl+ „. objects have m<]/m g ratio close to the value typical 
of normal PNe. If these nebulae were only partially ionized the 
true ntd/mg would be smaller than our estimate. Indeed, for H 1- 

11 Fig.l 1 3lis available online. 
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Fig. 6. The distribution of Sh/3 surface brightness for different types of GBPNe. Separate histograms have been constructed for 
objects with different ionization classes as described at the top of the panels. Above each bar the logarithm of the mean electron 
density of the PNe comprising it is given (in cm -3 ). 
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Fig. 7. Theoretical prediction of Sh/3 distributions for PNe with total gas mass of 0.2Mq, parameter e=0.75, and expansion at 
20km/s for the CSs of 0.57, 0.60, and O.64M0. Separate histograms are presented for objects with CS temperature T* <35kK (left 
panel), 35kK<T* <70kK (middle), and T* >70kK (right). The nebular electron densities are in addition plotted above the bar limits 
as expected from model calculations for PNe characterized by a given Sh/j value (log n e in cm -3 ). The numbers under CS mass give 
the cumulative percentage of objects expected in the analyzed temperature domain (in parenthesis if Sap is expected larger than the 
presented range of values). 



32 the ionized gas mass is exceptionally low m g =0.01Mo- For 
the other two bulge objects from the OC a+cr group, the derived 
m g is also only about one half of the typical mass of normal PNe. 

4.5. Dust temperature and infrared excess 

In Fig. [9] we present the dust color temperature derived from 
the 25 and 60jt/m IRAS flux ratios versus infrared excess IRE 
defined as 

IRE = F m I 223F m , 

where total infrared flux F/r has been estimated by computing 
the blackbody flux corresponding to the ratios of 25 to 60/mi 
fluxe s of IRAS. The values of bot h and IRE were adopted 
from Stasihska & Szczerba ( 1999). Figure |9] is found to be a 
good means of distinguishing all the different groups of GBPNe. 
There are two lines in this figure. The first represents the median 
IRE value for the normal GBPNe (dotted line) and the second 
the median T^ for normal PNe (dashed line). We note that al- 
most all the nebulae that we analyze here, DC cr , DC a+cr , and 
OC„ +cr , as well as PNe with emission-line CSs ([WR], WEL, 
and VL) are located above the latter line, meaning they all have 
dust hotter than the typical Tj of normal PNe. The possible 
reason for this is that they are presumably at a relatively early 



or intermediate stage of PNe evolution when the dust is hotter 
(Stasihska & Szczerballl999 l). In Fig. [6] we can see that the ma- 
jority of evolved objects (with log S///3<-2) are the normal PNe. 
In addition, selection effects may be important because strong 
dust emission (i.e., selection criterion for Spitzer/IRS targets) is 
more likely to be met in younger PNe with hotter dust when it is 
still relatively close to the CS. 

In Fig.|9j it is quite convincing that the OC a+cr belong to the 
objects with the hottest dust among bulge PNe. Apart from the 
above-mentioned relation with younger evolutionary age, other 
factors may be playing a role. The T^ can be higher in the 
case of low a md/m g ratio in the nebula and when amorphous 
silicate grains dominate the dust composition (see Fig. 7d,e of 
IStasiriska & S zczerball9 99). While there is only one OC a+cr ob- 
ject with definitely very low md/m g (M 2-23, see in Sect. 4.4) the 
members of OC n+cr are the only group of PNe analyzed here that 
show no clear signs of the C-based dust and only evidence of the 
silicate dust grains are present in their spectra. 

In their models, Stasihs ka & Szczerbal (1 19991) used the cir- 
cumstellar silicates dominated by amorphous grains, whereas in 
spectra of DC fr PNe and GBPNe with emission-line CSs there 
are clear signs of silicates in crystalline form only. It is not cer- 
tain how different forms of silicate grains may influence the 
dust temperature. In Fig. [9] one can note that both IC4776 - 
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Fig. 8. The mass of the dust versus the mass of the nebular 
gas for Galactic bulge PNe. The meaning of the symbols is the 
same as in Fig. [2] The horizontal short-dashed and vertical long- 
dashed lines indicates, respectively, median mj and m g for nor- 
mal bulge PNe (small dark grey symbols) and the dotted line 
indicate typical ma/nig ratio of these PNe. 



an unusual, unique WEL PN with 10/im amorphous feature and 
SwSt 1 - a [WR] type PN with a similar property belong to the 
hottest o f nebulae of their respective types (compare also with 
Fig. 3 of lGornv et aTll2001l) . On the other hand, the DC a+cr PNe 
that also exhibit amorphous silicates apparently have dust that is 
not convincingly hotter than for DC cr PNe that do not show the 
10/mi feature. 

The vertical dotted line in Fig.|9]indicates the median value 
of log IRE for normal PNe, which equals 0.63. This parameter 
and this particular value are important because it appears to sep- 
arate the different groups of GBPNe. First, to the right of this 
line are DC U+CT PNe (OC fl+c ,- are discussed below) and the pe- 
culiar VL PNe with cool CSs. In the case of IRE, the model 
calculations of IStasinska & Szczerbal d 19991) show that, unlike 
T f / discussed above, this parameter is less dependent on the evo- 
lutionary state of the PN. However, it depends strongly on the 
ntd/nig of the nebula as it is expected to be greater in objects with 
larger dust content. This is confirmed by the locations of DC fl+cr 
and VL PNe in Fig. [8] which indicate higher than average ra- 
tios of the dust-to-gas mass in these PNe. On the other hand, the 
DC cr objects characterized by m.d/m g at the normal level (just 
like [WR] PNe, or even like WEL PNe with md/m g below nor- 
mal) are located in Fig.|9]to the left of the dotted line with their 
IRE parameter being lower than average. 

By analyzing the IRE parameter, we find the most puzzling 
case to be the OC a+cr objects. The H 1-32 nebula with its rel- 
atively large md/m g is also characterized by a large IRE in ac- 
cordance with what we have just discussed. However, the other 
OC a+cr PNe with their (at most) moderate nid/m g should have 
much lower IRE values. This is not the case. The most extreme 
example is M 2-23, which has the lowest md/m g ratio yet its IRE 
value is among the largest. By analyzing the model calculations 



Fig. 9. The dust temperature versus the infrared excess param- 
eter for Galactic bulge PNe. The meaning of the symbols is the 
same as in Fig. [2] Dashed and dotted lines indicate median val- 
ues for normal bulge PNe (small black symbols). 

of Stasihs ka & Szczerbal (I 19991) . it seems that for OC I/+cr the fac- 
tor supporting larger IRE values may be their high density, e.g., 
due to thei r slower expansion. Accordin g to the models (Fig. 13a 
and 13c of Stasihsk a & Szczerball 999) IRE settles at an approx- 
imately constant level during an early phase of nebular evolution 
and becomes quite insensitive to other parameters. But if the 
nebular expansion is slow, as in the case of OC fl+c/ , the time be- 
fore it occurs is extended and IRE can remain large for a longer 
time. 

Finally, we should mention that the large IRE value of [WR]- 
type SwSt 1 nebula makes it again more similar to OC fl+CJ - or 
DC fl+cr than to bulge [WR] PNe in accordance with the large 
ma/trig ratio of this object that can be inferred from Fig. [8] 

4.6. Chemical composition 

We present the chemical composition of the analyzed PNe de- 
rived from optical spectra w ith the classi c al em pirical method. 
The values were taken from iGornv etaLl (|2009) where the ap- 
plied method have been described in detail. One import ant dif- 
ference of our study from analysis of Gorny et al. (2009) is that 
we do not limit our discussion of DC„ , DC fl+c ,-, and OCacr ob- 
jects to parameters with errors smaller than 0.3 dex as in the case 
of some PNe their spectra do not allow for that quality. 

The chemical elements in PNe are frequently divided into 
two groups. For one group of abundances it can be safely con- 
sidered that they remain mostly unchanged during the previous 
evolution of the CS and therefore the values found in PNe rep- 
resent the primordial abundances of the matter the progenitor 
star was born from. An example of such an element is oxygen, 
which is regarded as being mostly undisturbed since the object 
was born, at least in the case of PNe in metal-rich environm ents 
such as the Galactic bulge (see e.g.. IChiappini et al]|2009l for 
detailed discussion). The abundances of other elements are how- 
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Table 3. Abundances of oxygen and nitrogen adopted from calculations by iGornv et"a"D J2009) and compared to earlier 
published results. 





this work 


literature 




log U/rl+lZ 


lOg 1N/U 


(log U/rl + Iz, log IN/U) 


I 1 1- 1U 


O.OJ 


-1 OfS 

1 .UVJ 




Hb 6 


8 (\(\ 
o.w 


U. 1 o 


(8.71; -0.06)* 4 * 


M 2-10" 


o.u / 


-U.Z.O 


(8.78;-)° (8.86; -0.49)* (9.00; -0.55)* (8.90; -0.59) 5 ' 


H l-50~ 


8.66 


-0.37 


(8.65; -0.51) £ 


H 1-20 


8.61 


-0.10 


(8.97; -0.56) £ (8.45; -) G (8.97; -0.23)* (9.30; -0.67) Sfl (8.73; -0.44)* 


Th3-4 


8.22 


-0.11 




M3-38 


8.39 


0.19 


(8.37; 0.32)* (8.79; -0.61) Sn (8.54; -0.22)* (8.37; -) 5 ' 


M3-8 


8.38 


-0.26 


(8.59; -0.42) £ (9.02; -0.27)* (8.12; -0.27) c (8.15; -0.33) 5 ' 


H 1-40 


8.19 


-0.25 


(9.09; -0.04) £ (8.23; -0.45) HS (8.70; -0.62)* (8.62; 0.02f 


M2-23 


8.35 


-0.99 


(8.48; -1.44) £ (8.67; -1.00)" 4 ' (8.21; -0.13) 57 (8.22; -0.82)* (8.48; 0.20)* 


He 2-260 


8.00 


-0.76 




H 1-32 


8.59 


-1.21 


(8.55; -0.23) 5 ' (8.12; -0.53)* 


H 1-35 


8.65 


-1.19 


(8.31; -0.92)* (8.27; -0.89)* 



" The values ado pted are the average s of o riginal calculations (see Table 4) bylGornv et al. (2009) and from Stasiriska et al. ( 1998). 

Ref erences: AK -lAller & Kevesl d 19871), C -ICosta et alj dl996l). E -lExter et al J Bool). G -iGutenkunst et all (120081). R - iRatag et al J d!997l) 
K - iKoppen et a I dl99ll). Sa -ISamland et alJ d!992h. S t -IStasinska et alj dl99 
Note: Results bv lCuisinier et alj (2000). Escuder oetaT] d2004l) . lG6rnv et al.1 d20O 



^_HS -Ivan Hoof & van de Steend dl999h . 

and lWang & Liul (1200717 are not presented as the line measure- 



ments from these papers were used bv lGornv et al . (2009). 



ever expected to be changed as the result of various physical pro- 
cesses. Their abundance ratios are modified by nuclear reactions 
and mixing that can bring some freshly synthesized matter to the 
stellar surface in so-called dredge-up processes. An example of 
such an element is nitrogen. 

In Fig. [10] we plot the abundance ratios O/H versus N/O for 
the different types of GBPNe. These data can also be found in 
Table 3, where there are compared to earlier results from the lit- 
erature. As can be noted, the DC cr objects (except H 1-16) have 
locations in this plane that are compatible with the majority of 
bulge PNe, including normal ones and those with [WR] or WEL 
type CSs. The locations of DC a+cr are different as they show a 
substantial underabundance of oxygen related to both hydrogen 
and nitrogen. The OC a+cr PNe form another extreme with N/O 
showing the lowest values among GBPNe, while O/H seems nor- 
mal. Interestingly, the unusual objects IC4776 and SwSt 1 are 
located in the same region of the plot as OC a+cr PNe PI 

Before considering the possible reasons of this behavior, one 
must first discuss whether the obtained oxygen abundances are 
reliable. As directly indicated in Fig.QJJby the errorbars, the for- 
mal errors of derived abundances are sometimes very large, es- 
pecially for the DCa+c, PNe. They were calculated by propagat- 
ing the possible observational errors of measured spectral lines 
into comp uted parameters us ing the Monte Carlo method (for 
details see IGornv et al.1 12009). The typical errors for the other 
PNe presented in Fig.QJJare indicated with the errorbar cross in 
the bottom-right corner of the plot. 

Despite the large individual errors, we observe that the 
DC fl+cr PNe clearly have oxygen underabundances of more sim- 
ilar value than expected, because the random errors should pro- 
duce a greater scatter in the distributions than observed. In 
Fig. [10] we mark with a solid thick line the relation between 
O/H and N/O abundances for GBPNe with logN/O>-0.8 and 
good quality data (rejecting only some clear outliers). Assuming 
that this holds for all PNe and that observed deviations from it 
are caused only by individual errors the probability of finding 



all four DC f(+cr at their present locations in Fig.[10]can be evalu- 
ated as being smaller than 1%. This hypothesis can therefore be 
safely rejected meaning that the O/H vs. N/O relation for DC a+cr 
is truly significantly different from that of other GBPNe. Also 
using the Kolmogorov-Smirnoff 2D nonparametric test, the hy- 
pothesis that the difference between DC fl+c , and normal PNe in 
the Fig.[10]is statistically meaningful is confirmed at about 99% 
confidence level. That for the DC a+cr object with the smallest 
errors (M 3-38) we infer the same effect of oxygen underabun- 
dancq^ is aiso ver y important. 

In Cols. 3 to 9 of Table 4. are list the abundances of nitrogen 
and other elements for PNe investigated here and the median val- 
ues for normal, [WR], and WEL GBPNe. In Fig.QT] we present 
the histograms of the N/H ratio for the different types of GBPNe. 
Comparing the DC n+e) . PNe with normal PNe, one can see no en- 
hancement of nitrogen as the median values are almost the same. 
In the case of neon, sulphur, and possibly chlorine there seems to 
be no depletion of these elements relative to oxygen (see Fig. [12] 
available online). This means that whatever the reason for the 
lower oxygen abundance in DC a+cr PNe a similar decrease in 
the abundances of these three elements follows. Only in the case 
of argon is there possibly a difference at the level of O.ldex in 
the median abundance ratio of Ar/O between DC fl+e) . and nor- 
mal PNe. 

By analyzing the histograms of DC„. objects in Figs. [TTI and 
[T_2]there is no indication of any important difference in chemical 
abundances between this group and normal PNe. 

The situation of OC a+cr PNe is different. As can be deduced 
from the histogram in Fig.[TT] their underabundance of N/O is 
caused directly by the considerably lower number of nitrogen 
atoms. By analyzing the data for other elements in Fig. [12] we 
also find that the abundances of neon and argon seem lower than 
the average levels in normal PNe. On the other hand, no differ- 
ence is noticeable for the abundances of chlorine and sulphur. 



12 We computed chemical ab undances o f SwSt 1 and IC 4776 using 
line intensities from lPefia et alj d200ll) and lExter etaD d2004h 



13 In fact, for M 3-38 we have two independent, high quality spectra 
that provide very similar results. 
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Fig. 10. The nebular abundance ratios N/O versus O/H for the 
different groups of GBPNe with meaning of the symbols the 
same as in Fig. [2] The solid line marks the relation between O/H 
and N/O abundances (for PNe with logN/O>-0.8). The dashed 
vertical line marks median O/H of normal PNe. 
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Fig. 11. Distributions of N/H abundance ratio for the different 
groups of Galactic bulge PNe. The same notations apply as in 
Fig.0 



5. Discussion 

There has been a long accepted opinion that post-AGB objects 
can belong to only one of the two groups: either those charac- 
terized by an oxygen-rich environment or those surrounded by 
carbon-dominated matter. However, soon after the first obser- 



vations of Galactic PNe by ISO satellite were released it be- 
come clear that objects exist, that is [WR] PNe with late-type 
CSs, that simultaneously contain both carbon-based and oxygen- 
based dust. A number of hypotheses hav e been invoked to ex- 
plain this phenomenon (see the review in IPerea-Calderon et"al] 
2009). However, the dual-dust chemistry phenomenon remained 
to be viewed as a rare event among PNe that happens to only 
some objects. Furthermore, these objects were already acknowl- 
edged to be very peculiar because of the hydrogen-deficient 
composition of the CSs and strong WR- type stellar winds. 

The Spitzer observation s of GBPNe dGutenkunst et alj2008l: 
IPerea-Calderon eTa I] 120091) yielded a number of discoveries. 
The first surprise was a higher than expected percentage of 
PNe with dual-dust chemistry, which as we discussed in the 
Introduction mi ght have been partly c aused by the lower sen- 
sitivity of ISO. iGutenkunst et alj d2008l) argued that dual-dust 
chemistry may be related to the binary evolution and the exis- 
tence of circumbinary disks. They attributed the observed rate 
to a higher proportion of [WR] PNe in the Galactic b ulge, 
which has not been ac tually proven dGornv et al.l [2009), the 
IGutenkunst et alJ (|2008) sample containing only one such ob- 
ject. 

An alyzing their sample of GBPNe IPerea-Calderon et~a"D 
(2009) noted some other important results. The first was that 
dual-dust chemistry phenomenon is present not only in PNe sur- 
rounding late type [WC] CSs but in all [WR] PNe and also 
in all other PNe with emission-line CSs (WEL, VL) as well 
as ma ny PNe without emission-line CSs. IPerea-Calderon et"al] 
(2009) also found that the percentage of PNe with C-rich ma- 
terial is confusingly high in the Galactic bulge since there is a 
well known deficit of C-rich AGB stars in that environment. We 
note that all GBPNe that contain PAHs are dual-dust chemistry 
so urces. 

IPerea-Calderon et~ai] (120091) noted that the population of O- 
rich AGB stars in th e Galactic bulge, wit h the exception of ob- 
scured OH/IR stars ( Vanhollebekef l2007l) . do not show any in- 
dication of crystalline silicates in their spectra. Thus the high 
detection rate of dual-dust chemistry found in PNe cannot be ex- 
plained by long-lived O-rich (primordial or circumbinary) disks. 
That most of GBPNe cannot originate in the binary systems is 
indicated also by the results of Miszal ski et alj ([2009). They re- 
ported that among 300 analyzed GBPNe only 21 (about 12-21%) 
show the signatures of periodic variability in the OGLE-III data, 
which may be interpreted as being caused by close binarity. 

The low-mass O-rich AGB stars in the Galactic bulge can- 
not bring enough carbon into the envelope to produce C-rich 
AGB stars since in the high-metallicity environment the ef- 
ficiency of their third dred ge-up is very limited. Therefore, 
Perea-Calderon et al. (2009) proposed that the most plausible 
scenario for creating C-rich AGB stars seems to be a final ther- 
mal pulse on the AGB (or just after), which would produce an 
enhanced mass loss, capable of removing/mixing (sometimes 
completely) the remaining H-rich envelope and exposing the C- 
rich layers. It would also generate shocks responsible for silicate 
crystallization in the ejected circumstellar shell. 

For [WR] PNe there are many arguments that the change in 
the compos ition of the CS occurs p redominantly on the AGB or 
soon after dGornv & Tvl enda 2000). Since dual-dust chemistry 
in the Galactic disk PNe had been revealed in [WR]-type ob- 
jects it was natural to expect that a successful scenario should 
explain simultaneously both the unusual composition of JWR]- 
type CSs and the dual-dust chemistry of their nebulae (Gdrnv 
2008 ). The final thermal pulse at the end of the AGB suggested 
bv lPerea-Calderon e t al. ( 2009J) as a most plausible theory meets 
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this expectation. However, it is unclear why dual-dust chemistry 
is so widely observed among GBPNe and res t ricted to only some 
Galactic disk objects. iPerea-Calderon et all d2009l) pointed out 
that metallicity may be a possible explanation since the metal- 
licities of GBPNe are on average higher than those in the Milky 
Way disk. 

The properties of GBPN e with emiss i on-lin e CSs of various 
types were investigated by Gorn v et all d2009t) . In the present 
paper, we have therefore concentrated on those PNe that do not 
contain emission-line CSs. After inspecting the Spitzer spectra, 
we found that they can be divided into thred^l clearly separate 
groups DC cr , DC fl+c) , and OC a+CT . Objects of the first group, 
DC C) -, have simultaneously both carbon-based dust (PAHs) and 
oxygen-based dust (crystalline silicates). In the second group, 
DC a+cr , there are both PAHs and crystalline silicates, but in addi- 
tion there are also amorphous silicates. In the last group, OC fl+cr , 
there is only oxygen-rich dust but in both crystalline and amor- 
phous forms. Our primary discovery is that this classification 
based on the dust features in the infrared spectra is reflected in 
some other, more general properties of PNe. 

The DC„- PNe have infrared spectra that most closely resem- 
ble those of PNe with emission-line CSs. One can note the sim- 
ilar signatures of crystalline silicates at 23.5, 27.5, and 33.8/im 
and of PAHs at 6.2, "7.7", 8.6, and 11.3//m. At the same, time 
the chemical composition of DC cr nebulae cannot be distin- 
guished from that of the WEL objects and in general from the 
majority of GBPNe. In principle, one could consider whether 
DC„- are not an earlier evolutionary phase of WEL PNe when 
stellar emission lines are not yet visible but our analysis of their 
evolutionary state do not give strong support for such possibil- 
ity. It cannot however be completely ruled out that for some DC„. 
PNe, stellar emission-lines have not yet been discovered because 
of the quality of the spectra. 

For the possible links between DC tr and [WR] PNe, it has 
to be noted that the latter objects are definitely brighter (see 
Fig. 4) and have more massive nebulae (Fig. 6). Both differences 
could be explained if DC cr PNe are less evolutionary advanced. 
But in that case, it is hard to accept that relatively strong stel- 
lar emission-lines (expected spectral types should range from 
[WC11] to [WC7]) could remain unnoticed in the optical spec- 
tra. We can also add that large-scale turbule nt motions are char - 
acteristic of PNe around [WR]-type CSs dGesicki et al.1 120061) . 
however no information about them have so far been reported 
for any ofDC„. PNe. 

The objects belonging to DC a+cr and OC n+cr groups exhibit 
far more pronounced differences from the other GBPNe. The 
two groups are linked by showing evidence of amorphous sili- 
cates at lOfim and belonging to the densest PNe in the observed 
Galactic bulge population. However, the other observational re- 
sults imply that the evolutionary status of DC a+cr and OC a+cr 
must be completely different. 

The high densities of DC fl+cr PNe can be regarded as a sign 
of the relatively short time that has passed since they left the 
AGB but it is more likely caused by their nebulae being more 
clumpy. DC fl+cr are also characterized by the considerable excess 
in extinction that could be attributed to some source of internal 
extinction. Finally, the nebular gas has a peculiar chemical com- 
position with oxygen being underabundant relative to hydrogen, 
nitrogen, and possibly argon but preserving normal levels when 
compared to other elements. This seems difficult to interpret in 



14 There exist in the Galactic bulge also PNe with only signs of crys- 
talline silicates but as explained in Sect. 2 we did not have enough data 
to discuss them as a separate class here. 



the framework of the standard chemical evolution of PNe pro- 
genitors. If the O/H in DC a+cr PNe were to represent the pri- 
mordial oxygen abundance of the matter from which they were 
formed, it would favor older objects born before the interstellar 
matter was enriched in metals. In that case however, the CSs of 
DC n+c , PNe should be slowly evolving low mass objects and the 
N/O ratio should not be increased since the effective dredge-up 
of nitrogen occurs in higher mass AGB stars. In contrast, our 
analysis of the evolutionary status, indicate that DC a+cr have in- 
termediate and sometimes clearly higher mass CSs. 

The low O/H abundance accompanied by a higher than nor- 
mal N/O ratio is possible if the ON cycle of nuclear reactions 
was active in the progenitor stars. In that case the O/H ratio 
no longer reflects the primordial oxygen content. The ON cy- 
cle seems to work preferentially within more massive stars in 
low-metallicity environments. The examples are some PNe from 
LMC that show a clear anti-cor relation b etween O/H and N/O ra- 
tios dLeisv & Dennefeldl (120061) . see also lChiappini et al.1 (12009) 
and the discussion therein). On the other hand, the ON cycle 
should have no effect on other elements. For the DC a+a PNe 
discussed here, this is not the case as the abundances of neon, 
sulphur and chlorine closely follow the depletion of oxygen. 

Finally, in DC a+cr PNe some depletion of oxygen could be 
possible because it is being trapped in dust grains. The percent- 
age of oxygen removed in this way from the gas to the dust may 
be metallicity dependent and is limited by the amount of silicon 
available. More could be depleted by means of ice growth but 
this would give rise to recognizable features in the infrared spec- 
tra. The depletion may simultaneously concern not only oxygen 
but also some other elements such as sulphur. However, it should 
have no effect on noble gas such as neon. As the Ne/O ratios of 
DCa+cr nebulae are not enhanced, this hypothesis should also be 
ruled out. 

The OC fl+e ,- PNe are characterized by very small diameters 
and the highest densities but at the same time expand more 
slowly than other GBPNe. For this reason, their present evo- 
lutionary state indicates they should be related to intermediate 
mass CSs (»0.60Mq). However, this is inconsistent with the 
very low metallicity of their surrounding nebulae as seen in the 
low N/O ratio derived for these objects as well as of neon and 
argon. In contrast, low metallicity argues for the low mass CSs 
originating in lower mass progenitors that were created before 
the Galactic bulge was effectively enriched in metals^ 

The ev olution of a low mass C S can be accelerated by high 
mass-loss. Kudritz ki et alj ([1997) detected very strong winds 
from OC fl+cr object H 1-35 (the final numerical value has not 
been given). On the other hand, He 2-2 60 has normal, low m ass- 
loss at the level of O.45xlO~ 7 M /yr dHultzsch et alj|2007l) . In 
both cases the chemical composition of the CS is normal, i.e., 
they are not H-deficient. 

The possible solution to the puzzle of the relatively fast evo- 
lution in the nebular phase and the low metallicity of OC fl+cr 
could be their origin not from single stars but from a binary sys- 
tem including a lower mass star. This star would not normally be 
observed in the PNe phase because its post- AGB evolution is too 
slow to ionize the ejected gas before it disperses. However, as a 
result of the mass transfer from the companion at earlier evolu- 
tionary phases of the binary system (or even merging in extreme 
cases), the future CS could increase its mass and evolve more 
rapidly during its post-AGB phase. The PN would become visi- 



15 In that case, the O/H ratio may not repres ent the primordial com - 
position, as is the case of PNe in the LMC - see Chiappini et al. (2009) 
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ble, however, the nebular abundances of some elements possibly 
being characteristic of a lower mass (older) progenitor. 

The evolution of the AGB star in the binary syst em can have 
other consequences. As argued by de Marcc^ d2009t) for nebulae 
created from common-envelope binaries, the abundances of e.g. 
nitrogen and carbon should be statistically lower that in "nor- 
mal" PNe. This is because the interaction with the companion 
will cause the AGB star to terminate this phase of evolution ear- 
lier than of the single star. As a result, no effective dredge-up will 
take place associated with the thermal pulses at the tip of AGB. 
There will also be no final thermal puls e. The nebula M 2-23 was 
investigated by Miszals ki et al.l (|2009) but no photometric vari- 
ations were found. Nevertheless, the low N/O ratio is the char- 
acteristic feature of OC fl+cr PNe. We note also that the OC a+cr 
are the only GBPNe that clearly do not have C-rich dust. At the 
same time, not all the surrounding silicate grains have been crys- 
tallized. All these results imply that the formation of OC a+cr PNe 
differs from that of the rest of the Galactic bulge population. 

This leads us to the main conclusion of our work: there is 
clearly no unique road to the formation of PNe even in a uni- 
form environment such as the Galactic bulge. This is noteworthy 
since the stars in the Galactic bulge that we are now able to ob- 
serve in their PNe phase are expected to originate mostly from 
a single episode of star formation. Obviously, there are PNe in 
the Galactic bulge with very different CSs and different chemi- 
cal compositions. They are also characterized by different prop- 
erties of the dust as seen in DC er , DC a+cr , and OC fl+e) . groups 
investigated in this paper. Nonetheless, the simultaneous pres- 
ence of PAHs and crystalline silicates dominates in the GBPNe. 
Therefore, the scenario of final thermal pulse at the end of the 
AGB that changes both the stellar composition to C-rich and 
at the same time allows the crystallization of existing O-rich 
grains remains the most plausible possibility for the majority of 
GBPNe. However, it is not always effective or in different ways 
for different stars. 



6. Conclusions 

We have investigated PNe without emission-line central stars 
located towards the Galactic bulge that have peculiar infrared 
spectra acquired by Spitzer/IRS. Among these objects, we have 
found three separate groups divided according to their composi- 
tion of dust grains: 

- DC„- - dual-dust chemistry PNe with simultaneous presence 
of both carbon-based dust (PAHs) and oxygen-based dust 
(crystalline silicates); 

- DC fl+cr - dual-dust chemistry PNe with simultaneous exis- 
tence of PAHs and crystalline silicates as well as amorphous 
silicates; 

- OC a+cr - PNe characterized by oxygen dust chemistry with 
only oxygen-rich grains in both crystalline and amorphous 
forms. 

We have analyzed a wide range of different properties of these 
PNe. Our main results are: 

- We confirm that dual-dust chemistry is a common phe- 
nomenon of PNe in the Galactic bulge and can occur in ob- 
jects not related to emission-line central stars. 

- The Properties of DC„- PNe do not distinguish them clearly 
from the majority of other PNe in the Galactic bulge. They 
have intermediate or higher-mass central stars. Their infrared 



spectra closely resemble those of PNe with emission-line nu- 
clei. Some DC cr may be evolutionary related to the latter ob- 
jects or may have undiscovered emission-line central stars. 

- DC fl+cr objects belong to the densest PNe in the Galactic 
bulge. Their derived m c {/m g mass ratios and infrared excesses 
IRE are higher than average. There is a possibility of exten- 
sive internal extinction. DC fl+cr PNe have intermediate and 
higher-mass central stars. The chemical composition of neb- 
ular gas is peculiar as oxygen seems underabundant relative 
to hydrogen and nitrogen but not to other elements (except 
possibly argon). This composition of DC fl+cr PNe cannot be 
explained in the standard picture of AGB star chemical evo- 
lution. 

- The OC fl+ „ PNe are the only analyzed PNe not showing 
dual-dust chemistry. They have hottest dust temperature 
and highest infrared excess IRE. OC a + cr have also very small 
diameters and are among the densest PNe in the Galactic 
bulge. However, their expansion velocities are smaller than 
average and therefore their evolutionary status indicates that 
OC a+cr can have intermediate-mass central stars. In contrast, 
the surrounding nebulae show low metallicity with an un- 
derabundance of nitrogen, neon, and argon. The domination 
of oxygen-based dust indicates in addition a low abundance 
of carbon. We argue that their properties are in qualitative 
agreement with scenarios of PNe formation not from single 
AGB stars but from binary systems. 
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Tabl e 4. Chemic al abundances for analyzed PNe and median values for bulge [WR], WEL, and normal PNe samples (from 
iGornv et al .1120091) . 
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Fig. 12. Distributions of abundance ratios for different groups of Galactic bulge PNe. For the normal PNe, the median value along 
with the 25 and 75 percentiles are marked with three short vertical lines above the histogram. Total number of objects included are 
given below sample names. 
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Fig. 13. Theoretical prediction of Snp distributions for PNe of O.IMq and expanding at 12km/s in the case of central stars of 
0.57, 0.60, and O.64M0. The separate histograms are presented for objects with central star temperatures T* <35kK (left panel), 
35kK<T* <70kK (middle) and T* >70kK (right). 



